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МОРЕІ, ОЕ ІММОХОЗЕХМ5ОВ О5ІХС РІЕКЕКЕХСЕ КОСАТІОМ5 
ОХ ТНЕ НЕХАСОМХАЇ, І АТТІСЕ 


У роботі запропоновано модель імуносенсора, яка грунтується на системі різницевих рівнянь на гек- 
сагональній решітці. Уведено клас решітчастих різницевих рівнянь із затримками в часі для моделювання 
взаємодії антигенів-антитіл усередині імунопікселів. Побудова моделі грунтується на ряді біологічних при- 
пущень щодо взаємодії колоній антигенів та антитіл, а також дифузії антигенів. Для опису дискретних у 
просторі колоній, локалізованих у відповідних пікселях, використовується апарат різницевих диференціаль- 
них рівнянь на гексагональній решітці. 
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Вступ 

Стрімкий розвиток науки і техніки 
потребує появи нових методів детекції. 
Тому в науці та промисловості зростає ін- 
терес до біосенсорів. Біосенсори є альтер- 
нативою до загальновикористовуваних 
методів вимірювання, які характеризу- 
ються поганою вибірковістю, високою 
вартістю, поганою стійкістю, низьким 
відгуком і переважно можуть використо- 
вуватися лише високодосвідченим персо- 
налом. Біосенсори - це нова генерація 
сенсорів, що використовує у своїй конст- 
рукції біологічні матеріали, які надають 
високу вибірковість, селективність, точ- 
ність, дають змогу здійснювати швидкі 1 
прості вимірювання |1|. Біосенсори ха- 
рактеризуються високою ефективністю 1 
широко використовуються у харчовій 
промисловості |2), при захисті навколиш- 
нього середовища |3|, в оборонній про- 
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мисловості |4|, але найчастіше викорис- 
товуються у медицині |5-7| як інструмент 
для постановки діагнозів. У цілому сі- 
мейство біосенсорів ділиться на дві час- 
тини. Перша пов'язана з рівнем рецепто- 
ра до біологічного матеріалу, який вико- 
ристовується в його будові. Як рецептори 
можуть бути ензим, протеїн, порферін, 
антиген або антитіло. Друга частина біо- 
сенсорів обмежена до шару провідника, 
де біологічний ефект перетворюється на 
вимірювальний сигнал, який може бути 
електрохімічним, імпедансним, ампер- 
метричним, оптичним та ін. 

Постановка проблеми 

Математична модель імуносенсора 
повинна враховувати просторово-часові 
властивості пристрою, в якому викорис- 
товується детектор. Відносно просторо- 
вої організації досліджувана модель по- 
винна грунтуватися на певній дискретній 
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структурі, яка буде враховувати взаємо- 
дію пікселів імуносенсора. З точки зору 
часових змін, модель повинна описувати 
процеси, відомі як популяційна динаміка. 
Саме тому поставлено проблему розроб- 
ки математичної моделі імуносенсора, 
яка б враховувала цілий ряд біологічних 
припущень щодо основних компонент 
пристрою - колоній антигенів і антитіл, 
локалізованих у пікселях, дифузії колонії 
антигенів між пікселями та ін. 

Аналіз останніх досліджень і публікацій 

Останнім часом велика увага до- 
слідників привернута питанням розробки 
та використання сенсорів 1 біосенсорів (1- 
7). У роботі (5| наведено грунтовний ог- 
ляд теоретичних основ проектування біо- 
сенсорів. Прикладні аспекти використан- 
ня імуносенсорів обговорюються у (381). 

Важливе значення у функціонуванні 
біосенсорів має фізичне явище флуорис- 
ценції, про що йде мова в |3, 4, 91. У робо- 
ті П2| сформульовано основні задачі, 
пов'язані з дослідженням стійкості в біо- 
сенсорах. У роботах (1, 11) йдеться про 
проектування сенсорів, в основу роботи 
яких покладено перебіг хімічних реакцій. 
У роботі (11) для такого роду сенсорів бу- 
ло запропоновано математичне моделю- 
вання у класі решітчастих диференціаль- 
них рівнянь. У даній роботі такий підхід 
буде використано для імуносенсорів. 

У П3) викладено основні результа- 
ти щодо використання рівнянь популя- 
ційної динаміки. У роботі (10) розробле- 
но і досліджується спрощена модель 
імунної системи за допомогою диферен- 
ціальних рівнянь із запізненням, що буде 
використана в даній роботі. 

Мета дослідження 

Запропонувати модель імуносенсо- 
ра, яка грунтується на системі різницевих 
рівнянь на гексагональній решітці. 

1. Структура імуносенсорів та їх 
характеристики 

Серед великого сімейства біосенсо- 
рів імуносенсори є типовими сенсорами, 
що містять шар рецептора, який чутли- 
вий 1 селективний, включаючи імобілізо- 
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ваний біологічний елемент, наприклад, 
антитіло, антиген або хаптен, які є імуно- 
логічними рецепторами для молекул, які 
досліджуються. В імунсорі (імуносенсо- 
рі) відбувається реакція, яка грунтується 
на взаємодії між антитілом і антигеном 
або маленькими молекулами (хаптена- 
ми). Антитіла часто називаються імуно- 
глобулінами тому, що вони є протеїнами, 
які пов'язані з імунною системою. 

Імуноглобуліни використовуються 
імунною системою для ідентифікації та 
нейтралізації чужорідних об'єктів. Вони 
використовують властивості зв'язування 
антигенів. Антигени 1 антитіла можуть 
використовуватися у шарі рецептора в 
біосенсорах. Зменшення властивостей, 
які пов'язані з антитілами під час проце- 
су імубілізації антигену, використовуют- 
ься в конструкції шару рецептора, де 
антитіла відіграють функцію аналітів 
(молекул предметного детектування) |8). 

Молекули, які пов'язані з детекцією, 
забезпечують зв'язування антитіл з анти- 
генами, утворюючи складні конструкції. 
При цьому між антигенами і анти-тілами 
утворюються дуже сильні зв'язки з конс- 
тантою зв'язування Ка - 103" - 107". (9). 

2. Решітчасті диференціальні рівняння 

Решітчасті диференціальні рівняння 
ефективно використовуються в багатьох 
прикладних науках, таких як хімічні ре- 
акції, обробка зображень, матеріалознав- 
ство і біологія | 141. 

У моделях решітчастих диференці- 
альних рівнянь просторова структура має 
дискретний характер. Решітчаста динамі- 
ка широко використана в задачах (14-16|, 
оскільки середовища, у яких популяційні 
види живуть, можуть бути дискретними, 
а не неперервними. 

3. Решітчаста модель взаємодії 
антиген-антитіло для гексагонального 
масиву імунопікселів 

Для решітчастої моделі взаємодії 
антиген-антитіло для гексагонального ма- 
сиву імунопікселів використано матема- 
тичний опис за допомогою нелінійних 
різницевих рівнянь із запізненням. 
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Розглянемо модель імуносенсора на 
основі гексагональної решітки. При цьо- 
му для нумерації імунопікселів (1, ,К), 
і,),,ка- М,М, іч )-К 0 використову- 
ється кубічна система координат |17). 

Нехай У, ,, (г) - концентрація анти- 


генів, ЕК зо концентрація антитіл в 


імунопікселі (1,7,Ю); і,,,ко-М,У, 


іж)нкоад. 
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Модель грунтується на біологічних 
припущеннях для довільного імунопіксе- 
ля (ї,),К). 

Антигени детектуються, зв'язують- 
ся 1, нарешті, нейтралізуються антитілами 
з деякою ймовірнісною швидкістю у». 

Припускається, що коли колонії анти- 
тіл відсутні, колонії антигенів регулюються 
логістичним рівнянням із затримкою: 


У (21) зачв-бУ, к(п- т), , (п), (0) 


де 8 і б, - додатні числа, а г20 
означає затримку негативного відгуку ко- 
лоній антигенів. 

1. Вводиться константа народжуваності, 
В» 0 для популяції антигенів. 

2. Антигени нейтралізуються антитілами 
з деякою  Ймовірнісної швидкістю 
у»б. 

3. Популяція антигенів намагається до- 
сягнути деякої межі насичення зі 
швидкістю 6,20. 

4. Розглядаємо деяку дифузію антигенів 
з шести сусідніх пікселів (2--1, /,К-1), 
(і-1,)-1Ю, (,)-1,К-к1), 
(2-1, 7,1), (1-1,)-1Ю, 
(,)-К1,к-1), (рис. 1) зі швидкістю 


Піксель 


ву, Ю 


Піксель 
(1-1, , 91) 


дифузії РА", де р»д0 і Лл»б є від- 
станню між пікселями. 

5. Введемо деяку постійну смертності 
антитіл й, 20. 

6. У результаті імунної відповіді щіль- 
ність антитіл збільшується з імовір- 
нісною швидкістю ЛУ. 

7. Популяція антитіл прагне до деякого 
рівня насичення зі швидкістю б зо 

8. Імунна відповідь відбувається з дея- 
кою постійною затримкою у часі 
т»б. 


Піксель 
(ЧІ, ), Е-1) 


Рис. 1. Гексагональна решітка, яка пов'язує шість сусідніх пікселів у моделі 
імунопікселя з використанням кубічних координат 
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Таким чином, 
р р р 
1, 3, 5, 8, 9, 11 - и о; 2- нач) 4- ин У; 


р р р р 
6 - мин оУ, Де ман о; 10- иа о); 12 мн) 


На основі вищенаведеної інформа- 
ції запишемо математичну модель взає- 
модії антиген-антитіло із запізненням для 
гексагонального масиву імунопікселів, 
яка грунтується на добре відомій моделі 


Марчука (11-13) і використовує просто- 


ровий оператор 5, запропонований у 
17) (додаткова інформація с. 10). 


4У, зх 1) 2 
шо тд) «ва У ЧУ а) 
(2) 
4Б, зх 1) 
Гн з СИ т ПУ їж (ет)- б іа 9728мИ уві (г) 
Математична модель (2) задана по- 
чатковими функціями (3): 
У Фе 20, ЕЕ Фев (20, гє|- г), з) 
У 0 В (0)»0. 
Для гексагонального масиву вико- 
ристовується дискретна дифузія для 
просторового оператора: 
54 ї- РА" запою Ж Уа ї У ан ї У іа Є У ін ї У ін и лу, уд | (4) 
ік) з 
| і дЕє-М31,М-1, ічічкоб, 


Кожна колонія піддається впливу 
антигенів, вироблених у шести сусідніх 
колоніях, розділених рівними відстанями 

Використовуємо граничну умову 


У, -0 для пікселів (і,),К) таких, що 


і,,К 


іє |- У, М) м ) є|- М, М|м Ке|-М, М) та 
іжі)чКкоб. 


Визначимо фазовий простір С сис- 
теми (2)-(4), як Банаховий простір непе- 


рервних функцій ф:|-г,0)-з КУ з 
нормою Й - зир |/(г). 
є -т,0) 


Нехай 


С; АК рн СУ або, К-- М, Мі Кк 20) 


-У,М, 
іч | нКке0 


є невід'ємними і обмеженими на 
І-г,0), У, , . (0) 2 0,3, ,,(0)» 0). 


Виходячи з біологічного обгрунту- 
вання системи (2)-(4), бачимо, що почат- 
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кові умови розв'язків задовольняють 
ГАЛАН ру єС,. 


і, і,),к ЗАВННХ 

із у-нке0 
Можна показати, що функціонал у 
правій частині системи (2) є неперервним 


і задовольняє локальну умову Ліпшиця 
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щодо (У, зх (1 Б зх (1))), з косу. У Прос- 


іч )нкад 
торі С. Таким чином, виходячи з теорії 
функціонально-диференціальних рівнянь 
з кінцевим ен: П 8), для будь- 


якого || фо Урі здотуж, Є С 
іч )якед 
система (2) має єдиний розв'язок 
КУ, з Бі, з (2), зхобмомо» ЯКИЙ ЗАДО- 
іч яка 


вольняє початкову умову (3). Можна 
показати, що коли фе С, то розв'язок є 
позитивним, тобто 


У ФУ 20,5, 5 (1,фу 2 0, 
і, ),к -- М,М,і- )-Кко0 на інтервалі 
існування. 
Для дискретизації, дослідження 
стійкості та стабільності, використову- 


ється підхід, який розроблений у (19) для 
систем типу «хижак-жертва» у даному 
випадку антиген-антитіло. У роботі да- 


ний підхід розширений у випадку решіт- 
частої моделі з дифузією. 

У даній роботі використовуються 
наступні позначення: 


. символ і - т, п для деякого цілого 
і, т, п, тп означає 
і - т, т-і1,...п; 

. |е| означає найбільше ціле число 
функцій; 

. а" з5ир, ,-Я(п) та а ап о,» (п) 


для будь-якої обмеженої послідов- 
ності (а(п)); 

е 7 ї Кк) - безліч невід'ємних ці- 
лих чисел. 

Модель імуносенсора з викорис- 
танням різницевих рівнянь на гексаго- 
нальній решітці 

Система (2) без дифузії апроксиму- 
ється наступним диференціальним рівнян- 
ням з кусково-постійними аргументами. 


ау, 
о (роув Пи Кр ел прі) бум, (ел ел АФ, з (0, 
0) 
42, (г) 
Же -- Й -пуУ, би -Де/пр)-б ЗБ Ки иодЄ, (1) 
для іє Їп/в, (пчкт)п| пе/7". систему (5) по п, г), де і«(п ТТ), тоді 
Позначимо, що / пі- отримаємо: 
Іг/я|- ге 2". Проінтегруємо останню 
аУ, 
о - -(8- УБ 5, «ап -т)- б У, з «ар -т)ї (т), 
(6) 
емо -(- У, (плоті) б (п) 
РІ ЕП з ПАТ бе її (0) 
Позначивши в системі (6) 
У уп) ЗУ, з (п), вп) ЗВ з. (пЛ), 
отримаємо: 
У Фа У, у (пуехрів - у, (п- г)- б п г) (т) 
Е ее Е і (пуехрі- непу У «абтенручад, "Б Е «(Т) 
Враховуючи  2-»(п-1)й можна ретна аналогова неперервна в часі 


спростити систему (7), додавши дифузію 
до першого рівняння, отримується диск- 
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У апа) ЗУ, м(пехрів-уБ (пот) б гео у (п) 
Емі) З В у(прехріи, пу г) ду уп) 


Поведінка моделі імунопікселя з ви- 
користанням різницевих рівнянь на гек- 
сагональній решітці (8) може бути ін- 
шою, порівняно із системою решітчастих 
диференціальних рівнянь (2). Проблеми 
еквівалентності різницевих і диференці- 
альних рівнянь розглянуті в роботі |20). 
Еквівалентність траєкторій різницевих 1 
диференціальних рівнянь, які отримують- 
ся з допомогою прямого та зворотного 
перетворення Ейлера або центральних 
різницевих схем, може бути лише для не- 
великих часових відліків. Використовую- 
чи нестандартну схему, запропоновану 
Міккенсом |21|, встановлено динамічну 
узгодженість між дискретним часом 1 
неперервним часом моделі в |20|, 221. 

Результати чисельного моделювання 


Е йте- 1.2 


25 00 
іпдевх і 


Е йте-5.2 


іпдех | 


00 
іпдех і 


Е йте-9.2 
0 


25 00 50 
іпавх і 


(9) 


Розглядаємо модель (3) при М-4, 


В - дтіп |, знадкобНк му з кіп 7, 
пит 8 

ца, 5,205-Я, 

у п'є 

2 

ан зо па09, р-о2ту, 


Д «0,3 пт. 

Результат чисельного моделювання 
реалізований на рисунку 2, де зображені 
динамічні зміни в пікселях моделі імуно- 
сенсора з використанням різницевих рів- 
нянь на гексагональній решітці у випадку 
М 4 для популяції антигенів при зна- 
ченні запізнення в часі г-20. 


Е йте«-3.2 


маше 


00 
іпаевх і 


Е йтеє-7.2 


маше 
З 


-25 00 
іпавх ї 


Е йтес11.2 


маїце 
З 


00 
іпдех і 


Рис. 2. Динамічні зміни в пікселях гексагональної решітки у випадку 
М -4 для популяції антигенів при г -20 
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Аналізуючи динамічні зміни в пік- 
селях імуносенсора з використанням різ- 
ницевих рівнянь, слід зазначити, що при 
г «16 в імунопікселях зберігається стій- 
кий ендемічний стан, а при г2 16 спосте- 
рігаються періодичні зміни та перехід до 
нестійкого ендемічного стану. 

Висновки та перспективи подаль- 
ших досліджень 

У роботі запропонована модель іму- 
носенсора, яка грунтується на системі різ- 
ницевих диференціальних рівнянь на гек- 
сагональній решітці із запізненням. Побу- 
дова моделі грунтувалася на ряді біологіч- 
них припущень щодо взаємодії колоній ан- 
тигенів та антитіл, а також дифузії антиге- 
нів. Для опису дискретних у просторі коло- 
ній, локалізованих у відповідних пікселях, 
використано апарат різницевих решітчас- 
тих диференціальних рівнянь. 

Результати чисельного моделювання 
показали, що якісна поведінка системи сут- 
тєво залежить від часу імунної відповіді г. 

У подальших дослідженнях необ- 
хідно провести дослідження стійкості в 
моделі імуносенсора на основі різнице- 
вих диференціальних рівнянь на гексаго- 
нальній решітці із запізненням. 
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КЕ5ОМЕ 


А.5. 5уег5йиК 

Моде! ої Ше іпипипо5еп5ог оп (Пе 
Базі5 об іїНегепсе едпабоп5 оп а 
Ппехагопа! Пайісе 

І Фе уогк Ше пає! ої іттиповеп- 
5ог 15 ргоро5едй, уурісП 15 Ба5зед оп Ше 5у5- 
гет ої ФіНегепсе едиайопя оп а рехабопа! 
Іасе. Тре соп5ігисйоп ої Ше тодеі 15 
Баз5ед оп а пигабег ої Біоіобіса! аззитр- 
йоп5 абоці Ше іпіегасйоп ої соЇіопіе5 ої 
апйєепя апа апибоаіез, аз ме! а5 пе діНи- 
зіоп ої апійєепе. 

Матеїу, Феу аге Ше ГоПом/пе: 

- ме раме зоте соп5(апі Біг(Фга/е Гог апії- 
збеп роршайоп; апійєєпя аге пецігай7ед 
Бу апибоаіс5 муїВ зоте ргобабішу га/е; 

-  апйреп5 аге пецігай/ед Бу апибоаіе5 
уйіб 5оте ргобабійсу га/е; 

- оме Бауе 5опіе ФіНиш5іоп ої апібодіе5 
оп 5іх пеїєббогіпя ріхеїз5 ул дїйи- 
810п гаїе; 

- ме рауе 5оте соп5(апі дігібгаїе ої апіі- 
Боайіез; 

-  а5 а гезиіїс ої шатипе ге5роп5зе, ме Пауе 
іпсгеа5е ої деп5у ої апйбодіез мі 
зоте ргобабійу гаге; апабоду рориїіа- 
оп (епдз8 (о 50те сапуїпе сарасіїу мів 
зоте гаге; 

-  шптипе гезроп5е арреаг5 мії 5о0те 
сопуапі те 4еау. 

Тре птоз: ірогіапі ШФіпє 15 0 іаКе 
іпіо ассошпі 5зрацаПу Аїзсгеге сПагасіег ої 
Фе тоде!. Оп Ще оїег Бапа, у/е зе Ше 
ргедагог-ргеу тиодйеї їог Ше Цдезсгірйоп ої 
йпттипе гезроп5е їп еаср ріхе!. ТБе йте 
деіау ої іттипе ге5роп5е 15 декзсгібед Бу 
Фе аррагация ої деіау Фійегепсе едиацопе. 
Бог Фе Фезсгірйоп ої дїзсгеіе 5зрасез їп Фе 
зрасе ої фе соіопіез5, Іосайгед їп Ше 
соггезропдїпео ріхеіз, Пе аррагаїця ої 
ФіНегепсе едиайоп5 оп а Бехаєопа! Іакісе 
15 ц5ед. Рійегепсе едиацопя агі5е їп папу 
арріеад зибіесів, 5исп а5 спетіса! геаспоп, 
ітаєе ргосе55іпє, ппакегіа| 5сіепсе, ап 
ріоїогу. Ії Ше тодеіз ої іНегепсе 
едиайоп5, Ше 5рапа! 5(гисішге Баз а дї5сгеїе 
срагасіег, апа дФіНегепсе дЧупатіс5 Пауе 
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гесепйу Беєп ехіеп5іусіу ц5ед (о птіодеі 
ріоїодіса! ргобіетзя 5іпсе Ше епуіїгоптепі 
їп мПіср Ше 5ресіе5 роршіайоп Пуе5 тау Бе 
Фізстеїе Биї пої сопійпиоиц5. 5исПп арргоасп 
15 арргоргіаге Бог оиг ргобіет. 

Тре гезціїв ої питегіса! тлодейпя ої 
дупатіс срапее5 їп Ше ріхеї5 ої а піодсі ої 
ап іттпипобвеп5ог млі Ше ц5е ої дійегепсе 
едчайоп5 оп а Бехаєопа! Іайсе ої а 
роршайоп ої апієеп5 аге д1уеп іп Фе рарег. 

МУПпеп апаїугіпє Фе Фдупатіс срапее5 
іп Фе ріхеіз ої бе шпапипобвеп5ог ур Фе 
ц5е ої Фійегепсе едцайоп5, її уха5 сопсій- 
дед баг ії а 4еіау ої Іс855 Фап 16 їп 
йатипорірз, а 5іабіе епдаетіс 5іаіе 15 піаїп- 
каїпед, апа іп офФег са5е8, регіодіс спапее5 
апд а ігап5ійоп (о ап ип5кабіе епдетіс 5(а(е 
аге оббегуеа. 
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